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Experimental energy spectrum and B(E2) values in 68 Ni and 90 Zr indicate a double-magic char- 
acter in these neutron-rich nuclei with N or Z — 40. The data nevertheless do not show any 
pronounced irregularity in two-nucleon separation energy. To understand the underlying physics, 
we carry out both shell-model and mean-field calculations. The shell-model calculation can well 
reproduce all the observations. It is understood from the mean-field results for 68 Ni that the shell 
gap at N = 40 disappears due to dynamical correlations of the isovector J = pairing interaction. 
In 90 Zr, however, such a dynamic process with the J = pairing appears not important because of 
the strong contribution of the J > interaction. We study also level schemes in the Ni isotopes and 
N = 50 isotones. We predict a new band built on the state in both 68 Ni and 90 Zr. The states of 
this band are dominated by two-particle- two-hole excitations from the /p-shell to the intruder g g / 2 
orbit. 

PACS numbers: 21.10.Dr, 21.60.Cs, 21.60.Jz, 21.10.Re 



I. INTRODUCTION 

The study of nuclear shell effects away from the valley 
of stability is one of the current topics in nuclear structure 
physics. The most interesting aspects are how the well- 
known shell effects, such as the occurrence of magic num- 
bers and the shape-coexistence phenomenon 0, 0, 
manifest themselves in exotic mass regions where nuclei 
have unusual combinations of neutron and proton num- 
ber. There have been intensive discussions on the issue 
of weakening of shell effect in neutron-rich nuclei. For 
example, the spherical N = 20 shell gap for light nuclei 
disappears in neutron-rich isotopes, leading to strongly 
deformed ground states and large E2 transition proba- 
bilities between the 2\ state and ground state (Of ). By 
using the shell-model approach, it has been demonstrated 
H that the magic number at N — 20 vanishes due to the 
proton-neutron attraction between spin-orbit partners of 
maximum j. On the other hand, there have been sug- 
gestions p that the strong deformation effects around 
32 Mg are induced by dynamical correlations, such as the 
neutron pairing correlations. 

It has been found 0, S S by several experiments 
that the neutron-rich nucleus 68 Ni (Z = 28, N = 40) 
shows a double-magic character: a relatively large 2^ 
excitation energy and a small B(E2,0f — > 2+) value, 
which is comparable to the cases of double-magic nuclei 
16 0, 40 Ca, and 48 Ca. The nucleus 68 Ni lies far from the 
neutron drip line, and the neutron energy gap between 
the /p-shell and the g 9 / 2 intruder orbit ap pea rs to be 
sizeable at N — 40. It was discussed in Ref. (lOL flll| that 



in the early mean-field calculations, a distinct shell gap 
that exists in the N — 40 nucleus 68 Ni disappears when 
quadrupole correlations are taken into account. For 68 Ni, 
it is remarkable that this nucleus does not show a pro- 
nounced irregularity in two-neutron separation energy, 
as expected for a typical double-magic nucleus. It was 
suggested C2 that small B(E2,0^ — ► 2^~) value is not a 
strong evidence for the double-magic character. We may 
thus conclude that the double-magicity nature in 68 Ni is 
still controversial and remains an open question. 

In general, shell closure leads to spherical configura- 
tions for the ground state, while breaking of a magic 
shell can produce coexisting deformed states. An im- 
portant indication for the emergent deformation is the 
appearance of low-lying + bands. The deformed struc- 
ture occurs as a consequence of nuclear correlations, 
which excite nucleons from the closed shell to a higher 
shell. For exam ple, the typical double-magic nucleus 56 Ni 
(Z = N — 28) [ll Cilia is known to have two collec- 
tive bands with large deformations coexisting with the 
spherical ground band. Therefore, it is very interesting 
to examine theoretically whether such collective bands 
exist also in 68 Ni. 

Similar discussions would also apply to the neutron- 
rich nucleus 90 Zr, which has a closed Z = 40 proton sub- 
shell and a strong N = 50 neutron shell closure. This is 
an interesting case to study the persistence of the Z — 
40 stability. Recently, ener gy levels and B(E2) values 
in 90 Zr have been measured [l6j, which showed a double 
magic character: a relatively large 2^~ excitation energy 
and a small B(E2,0^ — > 2^~) value. However, this nu- 
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cleus does not indicate a pronounced irregularity in two- 
proton separation energy. Moreover, it is known that a 
low-lying 2 state exists at Z = 40 in the N = 50 iso- 
tonic chain. Hence we can expect to see excited bands in 
90 Zr but perhaps with different structure. 

In this paper, we study the magicity at N or Z = 40 
and structure of excited 2 bands in the neutron-rich 
nuclei 68 Ni and 90 Zr. To understand the physics in 
a systematical way, we perform spherical shell-model 
calculations for the Ni isotopes and N — 50 iso- 
tones. Conventional shell-model calculations in the 
(I/7/2, 2p 3 /2, 1/5/2, 2pi/ 2 , 139/2) shell space for N,Z = 
30 — 36 are not possible at present because of the huge 
dimension of configuration space; we need to restrict the 
model space to the 2p 3 / 2 , I/5/2, %Pi/ 2 , an d I39/2 orbitals 
(hereafter called the /pg-shell). Of course, neutron (pro- 
ton) excitations from the I/7/2 orbit to the fpg-she\\ can- 
not be neglected for 68 Ni ( 90 Zr) 0, E3- Nevertheless, 
after all we shall see that the variations in B(E2) in the 
nuclei around 68 Ni ( 90 Zr) can be understood in terms of 
valence neutrons (protons) in this restricted model space. 
For the Ni isotopes, we employ an effective interaction 
starting from a realistic neutron G-matrix interaction 
based on the Bonn-C NN potential (called VMS interac- 
tion) fvi\ . For the N = 50 isotones, we use two types of 
effective interactions: the proton part of the VMS inter- 
action and the effective interaction of Ji and Wildenthal 
(called JW interaction) |is|. 

The paper is arranged as follows. In Sections ITT1 and 
IIIII we present the numerical calculations and discuss the 
results for Ni isotopes and N = 50 isotones, respectively. 
Conclusions are drawn in Sec. IIVI 



II. NI ISOTOPES 

A. Magicity in 68 Ni 

Let us first review what the experiment has found for 
the Ni isotopes. In Fig. 1, the experimental B(E2, 0^ — > 
2+ ) value and the first excited 2 + energy are shown as 
a function of neutron number N for 58 ~ 70 Ni. With in- 
creasing N, B(E2) decreases quickly and becomes the 
smallest at N = 40. In contrast, changes in the 2+ en- 
ergy (£2+) are quite small for 58_66 Ni, but E 2 + jumps 
to a large value at N — 40. Thus, with a pronounced 
large E 2 + and a small B(E2, 0^ — > 2^), these data seem 

to suggest a subshell closure at N — 40 in 68 Ni. 

We now carry out shell-model calculations for the Ni 
isotopes. The shell-model Hamiltonian is written as 



H 



4 ^ 

ck/?7<5 



(1) 



where e a are single-particle energies and V a $^& two-body 
matrix elements. Since 56 Ni is taken as a core, the 
model space is restricted to the fpg-sheW for neutrons 
and protons are assumed to be inactive. The proton 
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FIG. 1: (Color online) Comparison between the calculated 
and experimental values of (a) B(E2, 0^ — > 2{ ) and (b) E 2 + 

for the Ni isotopes. The calculated values are denoted by 
open circles and the experimental data 0, 0, ITfll l2cjj by solid 
squares. 



core-excitations from 56 Ni are taken into account implic- 
itly by the effective two-body matrix elements and the 
proton contributions are estimated from the KB3 cal- 
culations |2l| in 48 Ca. The neutron effective charge is 
taken as e„ = 1.0 so as to reproduce the experimental 
B{E2,0+ -> 2+) of 68 Ni [13. We use the VMS interac- 
tion starting from a realistic neutron G-matrix interac- 
tion based on the Bonn-C NN potential. 

As one can see in Fig. 1, our calculations reproduce 
nicely the observed trends in B(E2,0^ — > 2±) and _E 9 + 

In particular, a large 2+ excitation energy and small 
B(E2, 0+ -> 2+) value at N = 40 are correctly obtained. 
It should be pointed out that the proton core excitations 
may significantly contribute to the excitation energy and 
to B(E2,0^ -> 2l) in 58 Ni, and thus it is difficult to 
absorb these effects into the effective interaction and the 
effective charges. In addition, the very recent observa- 
tion H3 indicates a large B(E2, 0+ -> 2+) value in 70 Ni, 
which exceeds the calculated one. Figure 2 shows two- 
neutron separation energy S^n and the difference between 
two-neutron separation energies 82m defined respectively 
by 

S 2n {Z,N) = B(Z,N)-B(Z,N-2), (2) 
8 2n {Z,N) = S 2n {Z,N)-S 2 n(Z,N + 2). (3) 
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FIG. 2: (Color online) (a) Two-neutron separation energies 
and (b) differences between the two-neutron separation en- 
ergies denned in Eq. (3). The exact shell-model results are 
denoted by open circles and the experimental data [IflLl22| by 
solid squares. 
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FIG. 3: (Color online) Spherical neutron shell structure. The 
HF single-particle energy levels e a predicted in the HF calcu- 
lations with the VMS interaction for the Ni isotopes. 



In Eq. 0, B(Z, N) is the binding energy taken as posi- 
tive values. The quantity S2n is known as the most sen- 
sitive and direct signature for a (sub)shell closure. Our 
shell model calculations reproduce well the experimental 
values of S2 n and S 2n - As can be seen in Fig. 2, S^n and 
#2n are smooth functions, and in particular, do not show 
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FIG. 4: (Color online) (a) Total energies and (b)correlation 
energies in the shell-model and the HF calculations with the 
VMS interaction for the Ni isotopes. Note that the absolute 
correlation energies become small around N = 40. 

any notable changes at N — 40. Thus, 68 Ni has a large 
E 2 + and a small B(E2, Of -> 2f ), but no irregularity in 
S2n and no strong peak in 8m ■ It is therefore very inter- 
esting to further look into the Sm and <5 2n results from 
the viewpoint of the magicity in 68 Ni. 

Let us analyze the shell-model results in Fig. 2 using 
mean-field procedures. We carry out Hartree-Fock (HF) 
and Hartree-Fock-Bogolyubov (HFB) calculations using 
the shell-model Hamiltonian (1). In the calculations, we 
impose spherical symmetry. The HF single-particle en- 
ergies are given by 



/3— occup 



(4) 



where $3fl=occup means the summation over the occupied 
states only. Figure 3 shows the HF single-particle ener- 
gies e a . The single-particle energy gap between g 9 / 2 and 
/p-shell varies from 4 MeV at N = 28 to 2.5 MeV at 
TV = 40, which shows a persistence of a large shell gap at 
this neutron number. As we shall discuss below, this gap 
in the static single-particle picture will be washed out by 
dynamic correlations. 

The total HF energy is expressed as 



E 



/3=occup 



(5) 
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On the other hand, the HFB approximation is carried out 
with the following procedure. The HFB transformation 
is given by 



u a c n 



(6) 



where a is the time reversed state to a and the occupa- 
tion numbers v 2 a satisfy the following equation 



1 



A 



(7) 



Here the self-consistent mean-fields, the self-consistent 
pairing gaps, and the canonical single-particle energies 
are respectively defined as 



and the total HFB energy [2^ is 
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HFB — 
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(9) 
(10) 
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The neutron chemical potential A is determined by the 
neutron number conservation 
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FIG. 5: (Color online) (a) Calculated two-neutron separation 
energies and (b) differences between the two-neutron separa- 
tion energies denned by Eq. (3) in the mean-field approxima- 
tion using the VMS interaction. The HF results are denoted 
by open triangles and the HFB ones by the solid squares. 
Note that at TV = 40 irregularity of S2n appears in the HF 
calculation and disappears in the HFB calculation. 



Eqs. and (|12f> are solved iteratively. In this paper, 
however, we get the solutions by minimizing the total 
HFB energy (|llf> with the neutron number conservation 
(|12f) under the normalization condition u\+ = 1. 

The total shell-model energies E$m and the HF ener- 
gies -Erf are plotted in Fig. 4(a), and the correlation 
energies, defined as E cor[ — Eqm — -Ehf, are shown in 
Fig. 4(b). The correlation energy exhibits a characteris- 
tic pattern where the absolute value is the largest at N — 
34 but has a local minimum at N = 40. The reduction 
in correlation energy at N = 40 would be attributed to 
the small pairing gap Ax/2 of the p\/i orbit with a small 
3- 

Calculations for two-neutron separation energy S% n 
and the difference between two-neutron separation en- 
ergies bi n are shown in Fig. 5. One can clearly see the 
irregularity in Si n and a peak in Si n in the HF calcu- 
lation for 68 Ni, which suggest a large energy gap and a 
subshell closure at N — 40. This result is consistent with 
the most of the Skyrme HF (SHF) and relativistic mean- 
field (RMF) calculations, which produced a distinct S^n 
peak at N = 40 pjj 1 1 lH . However, as seen in Fig. 5, 
the irregularity in Sin and peak in Si n do not show up 
in the HFB calculations when the T = 1, J = pairing 
interaction is included. We may therefore conclude that 



the T = 1 , J = pairing interaction is responsible for the 
observed smooth behavior in Sm and 8im and thus for 
the disappearance of a magicity character in 68 Ni. This 
conclusion is different from that of the SHF and RMF 
calculations in which the disappearance of the Sm peak 
is caused by quadrupole correlations 0, E| . 

The above conclusion is reinforced by the following 
analysis. To see the role of the T = 1,J = pairing 
interaction in the shell model calculations, we divide the 
two-body interaction Hi n t in the total Hamiltonian Q 
into two parts 



Hi, 



H J=0 + H J>0l 



(13) 



where Hj = q is the T = 1, J = pairing interaction and 
Hjyo = H — Hj-q. Figure 6 compares different calcula- 
tions for Sm and 5i n - We evaluate Si n and Si n by using 
the binding energy B(Z, N) calculated from the expec- 
tation values (H — Hj = q) and (H — -ffj>o), and compare 
them with the results of the full Hamiltonian. All these 
calculations use the same ground-state wavefunction ob- 
tained from diagonalization of the total Hamiltonian JQl . 
Now the significant role of the T = 1 , J = pairing in- 
teraction is clearly shown: when Hj = o is switched off, 
Sin exhibits irregularity and a large peak in Si n is seen 
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FIG. 6: (Color online) (a) Two-neutron separation energies 
and (b) differences between the two-neutron separation en- 
ergies defined by Eq. (3) from the shell-model calculations 
using the VMS interaction. The exact shell-model results 
are denoted by open circles, the expectation values neglect- 
ing the T = 1, J = interactions by open triangles, and the 
expectation values neglecting the J > interactions by open 
diamonds. Note that only the (H — Hj=o) result shows irreg- 
ularity at N = 40. 
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at N = 40, whereas in (H — iJj>o), no irregularity in S2n 
and no peak in #2n can be seen. 

In Fig. 7, we further examine the expectation values 
for various Hamiltonian terms. For the quantity (Hj = q), 
one sees that the contribution of the T = 1 , J = pairing 
causes a bending at N = 40. On the other hand, (Hj >0 ) 
increases monotonously with N. The total expectation 
value (Hi n t) m Fig. 7 corresponds to the correlation en- 
ergy in Fig. 4 (b). Thus we have understood the source 
of the seeming irregularity in 5^ and the peak in 52n 
(see Fig. 5). The irregularity shows up in two-neutron 
separation energy at N = 40 if the T = 1, J = pairing 
interaction is missing. Inclusion of the T = 1, J = pair- 
ing interaction washes out the irregularities in 5 2 n and 
&2n found in the HF calculations, and thus explains the 
observations. It was inferred from the discussion of the 
<7g/ 2 occupation number that the erosion of the TV = 40 
shell gap is attributed to the pairing correlations 0. 
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FIG. 7: (Color online) Expectation values of Hj = o and Hj > q 
in Eq. (13), which are denoted by open circles and open 
triangles, respectively. The total interaction energy (Hint) is 
also depicted by solid circles. Note that (Hj-q) displays a 
bending at N = 40. 



B. Level structure in 64 " 68 Ni 

In this section, we discuss the structure evolution along 
the isotopic chain 64_68 Ni. Figure 8 shows the experi- 
mental and theoretical level schemes. The B(E2) values 
have been measured only for the first transition between 
the 2+ state and the ground state [n|. Since in 64 Ni, the 
O2, 2^~, and 4^ states all lie around 2.7 MeV and their ex- 
citation energies are approximately twice the first excited 
l\ energy (~ 1.38 MeV), the level sequence appears to be 
consistent with that of an harmonic vibration. This se- 
quence is typical for low-lying excitations in spherical nu- 
clei. Anharmonicity of the 2-phonon states (0;j~, 2^~, 4+) 
becomes large in 66 Ni, and the harmonic pattern breaks 
down completely in 68 Ni where the 0^" level drops down, 
and appears below the 2+ level. 

We carry out shell-model calculations using the VMS 
interaction. The results are compared with data in Fig. 8 
and the predicted B(E2) values are summarized in Table 
I. The calculations can well reproduce the experimental 
energy levels and the B(E2, 0^ — > 2^) values, and the 
systematic behavior of the low-lying J state is also rea- 
sonably described. It is striking that in our results, an 
excited band is formed in 68 Ni based on the 0^~ state. 
The E2 transition probability B(E2,0^ -> 2^) is quite 
small in 64 ' 66 Ni, but becomes rather large in 68 Ni. The 
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values B(E2,A+ -> 2+) and B(E2,6f -> 4+) in 68 Ni are 
also large. In contrast, B(E2, Of — ► 2^~) in 68 Ni is found 
smaller than those in 64,66 Ni. Thus we have predicted a 
new band in 68 Ni, as shown in Fig. 8. 

In order to see how this band is formed, we compare 
neutron occupation numbers of the P3/2, /s/2>Pi/2) an d 
gg/2 orbits in the relevant low-lying states in 66 Ni and 
68 Ni. As one can see in Fig. 9, except for the 6f state, 
neutron occupation numbers in the low-lying states in 
66 Ni are dominated by the /p-shell components. This is 
because the Fermi energy of 66 Ni lies below the pi/% orbit. 
However, neutron occupation numbers in 68 Ni show very 
different values ||. Except for the ground state, occu- 
pation number of the gg « orbit in all low- lying states in 
68 Ni increases by more than two units. This means that 
two neutrons are excited from the /p-shell to the g^ii 



TABLE I: B(E2) values for the positive-parity yrast states 
and some excited states in 66 Ni and 68 Ni. The calculated 
values are the shell-model results using the VMS interaction. 
Data are taken from Refs. [fl l9Hl9|. 
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FIG. 9: (Color online) Neutron occupation numbers of the 
/pg-shell orbits for the low-lying levels in (a) 66 Ni and (b) 
68 Ni. 



orbit in these states in 68 Ni [Til l24j. To see the struc- 
ture of the low-lying states more clearly, we calculate the 
probability of n-particle-n-hole (np-nh) excitations from 
the /p-shell to the gg/2 orbit, defined by 



Pn = 



(Nn) 



(14) 



where N n are the np-nh operators from the /p-shell to 
the orbit. Table II lists the probabilities of np-nh 
excitations in the relevant low-lying states of 66 Ni and 
68 Ni. In 66 Ni, in all low-lying states except the state, 
the dominant components are the Op-Oh excitations but 
with considerable mixing of the 2p-2h excitations. The 
4p-4h excitations are quite small in these states. The 
6^ state in 66 Ni has almost a pure 2p-2h component. In 
contrast, the low- lying excited states in 68 Ni show very 
different structures. While the ground state has mixed 
2p-2h and Op-Oh components with nearly equal probabil- 
ity, the low-lying excited states have mainly the 2p-2h 
component with considerable mixing with the 4p-4h ex- 
citation. A large E 2 + and a small B(E2, 0^ — + 2+) in 

68 Ni (see Fig. 1) would be alternatively explained as fol- 
lows. Once the odd-parity fp orbits are filled at N = 40, 
at least two- neutrons have to jump to the intruder g 9 / 2 
orbit to create a 2^~ state, and therefore the energy E 2 + 

increases 0,0,13- The E2 transition between 2+ and 
0± in 68 Ni becomes small just because the two states have 
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TABLE II: Probabilities of n particle-hole excitations for the 
low-lying states of 66 Ni and 68 Ni. 
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different structure. Interestingly, we indeed see from our 
calculation that the band is built on the 2 state. This 
happens because all the excited states belonging to this 
band have a similar structure with the 2p-2h excitations. 



In order to visualize the shape of 68 Ni, we use the CHF 
method with the following quadratic constraint flh\ 

H' = H + a]T((Q 2M }- gM ) 2 + /?(<J,)-^) 2 ,(15) 

where Q 2p _ and J x are the isoscalar quadrupole oper- 
ators and the ^-component of angular momentum op- 
erator, respectively. The q p s are constant parameters: 

1o = \flh ( l cos ~f, q ±2 = y^f^qsin7, and q ±1 =0, where 
q is the isoscalar intrinsic quadrupole moment and 7 is 
the triaxial angle. We set j x = \J J {J + 1) with J the 
total angular momentum of the state. The parameters, a 
and /3, are taken so as to achieve a convergence for an it- 
eration calculation with the gradient method. Then, po- 
tential energy surface (PES) is defined as the expectation 
value {H} with respect to the CHF state for given q and 
7. Figure 10 shows the contour plot of the PES in the (7-7 
plane for 68 Ni. We find that the PES minimum exhibits a 
spherical shape and an oblate softness. This is consistent 
with our previous discussions on the shell- model results, 
namely, a large E 2 + and a small B(E2 1 Of — > 2+) in 68 Ni. 

The PES figure in Fig. 10 is in contrast to the charac- 
teristic feature of an oblate-prolate shape-coexistence in 
68 Se US- 

III. TV = 50 ISOTONES 
A. Magicity in 90 Zr 

In the previous section, we have discussed several un- 
usual properties found in 68 Ni, which are associated with 
the subshell closure at N — 40. A related question is 
how neutron-rich nuclei with Z = 40 behave. Figure 11 
shows the experimental B(E2, Of — > 2f ) and the first 
excited 2f energy as a function of proton number Z, 
for some N = 50 isotones. For both B(E2, Of -> 2f) 
and E 2 + values in Fig. 11, we find remarkable similar- 
ities as seen in Fig. 1: with increasing proton number 
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FIG. 10: Contour plot of PES on q - 7 plane in the CHF 
calculation for 68 Ni. 

Z, B(E2) quickly increases until Z = 34 and then de- 
creases from Z — 36 to Z = 40. The first excited 2f 
energy E 2 + goes up gradually and is peaked at Z = 

40. Again, in terms of B(E2,0+ -> 2+) and E 2 +, 90 Zr 
seems to be a double-magic nucleus. It should be pointed 
out that B(E2, 4 t — > 2+) shows different behavior from 
B(E2,0f — > 2^") |2Sf. Recent lifetime measurements for 
96 Pd and 94 Ru corroborate the tendency of this behavior 
for N = 50 [23| • Moreover, it has been shown recently 
that the exact strengths for these transitions cannot be 
reproduced in a T = 1 model space but require neutron 
excitations across the N = 50 shell |30j . 

We carry out shell-model calculations for the N = 50 
isotones. Since 78 Ni is taken as a core, the model space 
for proton is restricted to the /pg-shell, and the neutrons 
are assumed to be inactive. The proton effective charge 
is taken as e p — 1.8 for the VMS interaction and e p = 2.0 
for the JW interaction so as to reproduce the experimen- 
tal B(E2, Of -> 2+) value of 90 Zr 17}. We use two types 
of effective interactions: the proton part of the VMS in- 
teraction and the JW interaction. As one can see in Fig. 
11, the calculations nicely reproduce the observed trends 
for both B(E2, Of -> 2+) and E 2 +. 

Figure 12 shows the two-proton separation energy S^p 
and the difference between two-proton separation ener- 
gies 82p for this isotonic chain, defined by 

S 2p {Z,N) = B(Z,N)-B(Z-2,N), (16) 
S 2p {Z,N) = S 2p (Z 1 N)-S 2p (Z + 2,N). (17) 

The experimental data do not show a signature for a sub- 
shell closure in 90 Zr, since no irregularity in S 2p can be 
seen. The shell-model calculations reproduce well the ex- 
perimental values of S 2p . In particular, the small peak 
in S 2p at N — 38 is well described. To understand these 
results, we analyze the role of the T = 1, J = pairing 
interaction (Hj—o) and the other interactions (-Hj>o) in 
the Hamiltonian, as done in the previous section (see 
Eq. I|13[)). In contrast to the case of the Ni isotopes, 
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FIG. 11: (Color online) Comparison between the calculated 
and experimental values of (a) the B(E2,0^ — » 2+) and (b) 
E 2 + for N = 50 isotones. Data are taken from Refs. [Trl YWi 

|23|. The shell-model calculations are carried out using the 
VMS and JW interactions. The calculated results are denoted 
by open circles (triangles) for the VMS (JW) interaction, and 
the experimental ones by solid squares. 
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FIG. 12: (Color online) Two-proton separation energies in 
the shell-model calculations with the VMS interaction. The 
exact shell model results are denoted by open circles, the ex- 
pectation values neglecting the T = 1, J = interactions by 
the open triangles, and the expectation values neglecting the 
J > interactions by the open diamonds. Experimental data 
fl9ll22| are denoted by solid squares. 
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Figs. I12f a) and Fig. I13f a) indicate that, while the 
T = 1, J = pairing interaction scarcely contributes 
to S2 P , the remaining interactions -ffj>o increases S^p 
significantly. Thus, the J > interactions are more im- 
portant for two-proton separation energy in the N = 50 
isotones. The -ffj>o contribution, however, does not pro- 
duce any notable irregularity in S2 P - For 52 P , we can see 
some differences between the VMS and JW interactions 
in Figs. I12f b) and ll^f b). Moreover, the Hj = q and -ffj>o 
contributions to S2 P in the VMS results are different from 
those in the JW results. 

FigurelTHshows the expectation values of Hj—q, -ffj>o 
and the total interaction energy (Hi nt ). Comparing Fig. 
rm with Fig. [7J we find that in the N = 50 isotones, 
(i?j>o) increases drastically with increasing proton num- 
ber, and becomes dominant when Z is large. There is 
no clear bending at Z = 40 in either curve (Hj=q) and 
(i?j>o). Thus these detailed results have explained the 
trends of two-proton separation energy in Fig. I12f a) and 
Fig. Ula). 

Let us now study the contributions from the above 
interactions to HF single-particle energies e a in the HF, 
HF+BCS, and HFB treatments. We also evaluate two- 
proton separation energy S2 P within these treatments. 



FiguresirSlandirplshow respectively the results calculated 
with the VMS and JW interactions. It is seen that in 
the VMS results shown in Fig. I15f a). the single-particle 
energy gap between gg/% andpi/2 orbits decreases quickly 
with increasing proton number. This causes a smooth 
variation in S2 P as seen in Fig. I15f b). All the HF type 
calculations do not produce irregularity in S^p. In the 
JW results in Fig. iTpT a). the single-particle energy gap 
between the gg/2 and pi/2 orbits remains large up to Z = 
36, but becomes small after Z = 38. The Fermi energy, 
as found in all HF type calculations, lies in the /p-shell 
for Z = 30 — 38, and between g$/2 & n d P1/2 for Z = 40 — 
46. Therefore, protons do not encounter a large energy 
gap when they are excited. Therefore, also with the JW 
interaction, irregularity in S^p is not produced (see Fig. 
I16f b)). We note that in both Figs. H5land lT^l the proton 
separation energies in the HFB calculations deviate from 
those of the other calculations when Z is large. Similar 
trend is obtained in the shell-model calculation without 
the J > interaction. 
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FIG. 13: (Color online) Same as Fig. 1121 except that the 
calculations are performed by using the JW interaction. 

B. Level structure in S6 Kr, 88 Sr, and 90 Zr 

Experimental level schemes for 86 Kr, 88 Sr, and 90 Zr 
are shown in Figs. 1171 andH%l B(E2) in these nuclei has 
been measured only for transitions between the 2^~ state 
and the ground state. Among the three isotones, the 0^~ 
level in 90 Zr is the lowest in energy and lies below the 
l\ state. We perform shell-model calculations using the 
VMS and JW interactions, and the results are compared 
with data in Figs. 1171 and Upland the B(E2) values are 
summarized in Tables III and IV. 

The calculations can reproduce the experimental en- 
ergy levels and B(E2, 0+ — * 2+) values. In particular, the 
systematical behavior of the 0^" state is well described. It 
is striking that the results show again an excited band in 
state. For the E2 transition proba- 
) , both calculations indicate a quite 
small value in 8B Kr and 88 Sr, but a very large one in 
90 Zr. Moreover, B(E2,4+ -> 2+) and B(E2,6f -> 4+) 
are found large in 90 Zr. In contrast, B(E2, Of — > 2f ) in 
90 Zr is smaller than that in 86 Kr and 88 Sr. All of these 
suggest strongly a new band in the Z = 40 nucleus 90 Zr. 



90 Zr based on the 0^" 
bility B(E2,0+ -> 2- 




28 30 32 34 36 38 40 42 44 46 48 



FIG. 14: (Color online) Expectation values of Hj = o and Hj > r ) 
defined by Eq. (13) for (a) the VMS interaction and (b) the 
JW interaction. (Hj = o) and (Hj>o) are denoted by open 
circles and open triangles, respectively. The total of them 
({Hi n t}) is also depicted by solid circles. Note that the con- 
tributions of (Hjyo) becomes large with increasing proton 
number. 



in two shell-model calculations with the VMS and JW 
interactions. The results for 88 Sr and 90 Zr are listed in 
Tables V and VI, respectively. For 88 Sr with the VMS 
interaction (Table V), the ground state and the 2f state 
have a dominant component of the Op-Oh excitation, and 
the f and 4^ states have comparable probabilities of the 
Op-Oh and 2p-2h excitations. Note that the 2f state in 



TABLE III: B(E2) values for the positive-parity yrast states 
and some excited states in 88 Sr and 90 Zr. Data are taken 
fromRefs. |T1.IT^. The calculated values are the shell-model 
results using the VMS interaction. 



To confirm the above findings, we further study the 
probability of the np-nh excitations defined by Eq. (14), 
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FIG. 15: (Color online) Spherical proton shell structure, (a) 
The HF single-particle energies €-a obtained by the HF cal- 
culations with the VMS interaction and (b) the two-proton 
separation energies using the mean-field calculations for the 
N = 50 isotones. The HF results are denoted by open trian- 
gles and the HFB ones the solid squares. 



88 Sr can be made by lp-lh excitations from (/5/2,J>3/2) 1° 
P1/2, which contribute to the E2 transitions. In the JW 
results (Table VI), the 2p-2h components are dominant 
in the 0^ and 4+ states. The 6^ state is almost a pure 2p- 
2h excitation in both VMS and JW interactions. One can 
thus expect that only the E2 transition B(E2, 0^ 2~\) 
is enhanced in 88 Sr. In 90 Zr, on the other hand, the 
ground state has the Op-Oh and 2p-2h components with 
nearly equal probability and the dominant components 
in the 2^, 0^, and states are the 2p-2h excitation 



TABLE IV: Same as Table III, except that the calculations 
are performed by using the JW interaction. 
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FIG. 16: (Color online) Same as Fig. 1151 except that the 
calculations are performed by using the JW interaction. 



TABLE V: Probabilities of np-nh excitations in the low-lying 
states for 88 Sr and 90 Zr, where the VMS interaction is used 
in the shell-model calculations. 
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0.142 



mixed with the Op-Oh component. Similar results are 
found in both calculations. From this analysis, we can 
understand that an excited band is formed on the 0^ 
state in 90 Zr because the 0^", 2^~, 4^", and states all 
have a similar structure with a large component of the 
2p-2h excitations. 



IV. CONCLUSIONS 

We have studied the magicity of TV or Z = 40 and the 
level schemes for the neutron-rich nuclei 68 Ni and 90 Zr 
by means of the shell-model and the mean-field approxi- 
mations. For both nuclei with either N = 40 or Z = 40, 
their two-nucleon separation energies do not show any ir- 
regularity along the respective isotopic or isotonic chain, 
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FIG. 17: (Color online) Comparison between the experimen- 
tal and calculated level scheme for 86 Kr, S8 Sr, and 90 Zr. Data 
are taken from Refs. [l6l Il9||. The shell model calculations 
are carried out by using the VMS interaction. The widths of 
the arrows denote relative values of B(E2). The numbers by 
the arrows are the B(E2) values in Weisskopf units. 
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FIG. 18: (Color online) Same as Fig. 17, except that the 
theoretical results are obtained by using the JW interaction. 



in spite of the apparent double-magic feature shown with 
a comparatively large 2^~ excitation energy and a small 
B(E2,0f — > 2+) value. The reason why the separation 
energy does not exhibit irregularity has been found dif- 
ferent for the Ni isotopes and the N — 50 isotones. From 
the shell-model calculations using the VMS and JW in- 
teractions, we have suggested that the T = 1, J = 
pairing interaction is responsible for the absence of any 
irregularity in separation energy in 68 Ni. The irregularity 
appears in the HF treatment but disappears in the HFB 
treatment. This indicates that the shell gap at N = 40 



TABLE VI: Same as Table V, except that the calculations are 
performed by using the JW interaction. 
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disappears due to dynamical correlations of the isovector 
J = pairing interaction. In the case of 90 Zr, however, 
irregularity in two-proton separation energy does not ap- 
pear in the HF calculations. For the N = 50 isotopes, 
the J > interactions contribute significantly to the two- 
proton separation energy. 

We have also studied level schemes for 68 Ni and 90 Zr. 
We have predicted an excited band built on the 2 state 
in both nuclei. The dominant component of this band 
has been determined as the 2p-2h excitations from the 
/p-shell to the intruder g 9 / 2 orbit. The structure of the 
excited states of this band is quite different from that 
of the ground state. This happens because the opposite 
signs of parity between the gg/ 2 orbit and the /p-shell 
do not allow lp-lh excitations [24]]. The first excited 2f 
state in 68 Ni and 90 Zr lies higher, and B(E2, 0+ -> 2+) 
is relatively weak. The difference in parity between the 
fp and the gg/ 2 subshells leads to a small probability of 
quadrupole excitations across N = 40, and the large en- 
ergy gain due to pairing correlations in the g g / 2 subshell 
is responsible for the high 2 + energy in 68 Ni. 
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